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ABSTRACT: During embryo implantation in pigs, the uterine endometrium undergoes dramatic morphological and functional 
changes accompanied with dynamic gene expression. Since the greatest amount of embryonic losses occur during this period, it is 
essential to understand the expression and function of genes in the uterine endometrium. Although many reports have studied gene 
expression in the uterine endometrium during the estrous cycle and pregnancy, the pattern of global gene expression in the uterine 
endometrium in response to the presence of a conceptus (embryo/fetus and associated extraembryonic membranes) has not been 
completely determined. To better understand the expression of pregnancy-specific genes in the endometrium during the implantation 
period, we analyzed global gene expression in the endometrium on day (D) 12 and D15 of pregnancy and the estrous cycle using a 
microarray technique in order to identify differentially expressed endometrial genes between D12 of pregnancy and D12 of the estrous 
cycle and between D15 of pregnancy and D15 of the estrous cycle. Results showed that the global pattern of gene expression varied with 
pregnancy status. Among 23,937 genes analyzed, 99 and 213 up-regulated genes and 92 and 231 down-regulated genes were identified 
as differentially expressed genes (DEGs) in the uterine endometrium on D12 and D15 of pregnancy compared to D12 and D15 of the 
estrous cycle, respectively. Functional annotation clustering analysis showed that those DEGs included genes involved in immunity, 
steroidogenesis, cell-to-cell interaction, and tissue remodeling. These findings suggest that the implantation process regulates differential 
endometrial gene expression to support the establishment of pregnancy in pigs. Further analysis of the genes identified in this study will 
provide insight into the cellular and molecular bases of the implantation process in pigs. (Key Words: Pig, Uterus, Endometrium, 
Implantation, Microarray) 



INTRODUCTION 

Embryonic loss reaches up to approximately 20% by 
day (D) 20 of pregnancy in pigs (Bennett and Leymaster, 
1989). The most critical period in pigs is considered to be 
D12 of pregnancy, the time when implantation of the 
conceptus (embryo/fetus and associated extraembryonic 
membranes) begins and matemal recognition of pregnancy 
occurs via estrogen secretion by the conceptus (Bazer and 
Thatcher, 1977; Geisert et al., 1982). It has been suggested 
that asynchrony of embryo development and inappropriate 
interactions between the conceptus and the uterine 
endometrium are major reasons for embryonic mortality 
during early pregnancy (Geisert and Yelich, 1997). Thus, to 
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increase the rate of successful pregnancies, it is essential to 
imderstand the mechanism of the estabUshment of 
pregnancy. 

The implantation process is a well-coordinated 
interaction between the developing conceptus and the 
maternal uterus (Geisert and Yelich, 1997). At the time of 
implantation, the conceptus produces estrogen, which acts 
as a signal for maternal recognition of pregnancy in pigs 
(Bazer et al., 1986). It has been suggested that estrogen 
reorients prostaglandin (PGF2a) secretion from an 
endocrine to an exocrine manner and exerts an anti- 
luteolytic action (Spencer and Bazer, 2004). In addition to 
estrogen, the conceptus produces cytokines, such as 
interleukin (IL)-ip, interferon (IFN)-y and IFN-5, growth 
factors, and proteases (Tuo et al., 1996; Geisert and Yelich, 
1997; Lefevre et al., 1998; Ross et al., 2003a). In response 
to these factors, the uterine endometrium undergoes 
dramatic morphological and functional changes to become 
receptive to the developing conceptus. 
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In the presence of a conceptus dmng early pregnancy, 
the expressions of several genes in the uterine endometrium 
are known to change. For example, fibroblast growth factor 
7 (FGFT) is abundantly expressed in uterine endometrial 
epithelial cells during early pregnancy and is increased by 
estrogen treatment in endometrial explant cultures (Ka et al., 
2001). Secreted phosphoprotein 1 (SPPl, or osteopontin), 
which binds to integrin receptors and functions in cell-to- 
cell attachment, is localized to endometrial luminal 
epithelial cells during early pregnancy and is induced by 
estrogen (Garlow et al., 2002; White et al., 2005). In 
addition, lysophosphalidic acid receptor 3 (LPAR3), which 
is a receptor for lysophosphatidic acid that possesses a 
growth factor-like activity, is expressed in uterine 
endometrial epitheUal cells during early pregnancy and is 
induced by estrogen (Seo et al., 2008). Although it has been 
studied on the expression of many genes in the uterine 
endometrium during the estrous cycle and pregnancy, the 
pattern of global gene expression in the uterine 
endometrium in response to the presence of the conceptus 
has not been completely determined in pigs. 

There have been many efforts to understand the 
expression and function of the uterine endometrial genes 
responsible for the estabUshment and maintenance of 
pregnancy using various gene expression analysis 
techniques in pigs. Suppression subtractive hybridization 
(SHH) has been used to determine novel gene expressions 
in Meishan-Landrace conceptuses and endometrial tissues 
at D15 of pregnancy (Vallee et al., 2003). In a recent study, 
we characterized several differentially expressed genes in 
the uterine endometrium on D12 of the estrous cycle and 
pregnancy using an annealing control primer (ACP)-based 
PGR (AGP-PGR) technique (Ka et al., 2009). Microarray 
technology has also been used to analyze gene expression in 
uterine endometrial tissues on D30 (Ka et al., 2008) or in 
cyclic uterine endometrium in pigs (Lee et al., 2005; 
Whitworth et al., 2005; Green et al., 2006). However, there 
are Umited reports that have measured the global expression 
of genes in the uterine endometrium during embryo 
implantation in pigs. 

Therefore, this study aimed to determine the global 
expression of uterine genes at the time of embryo 
implantation using a microarray technique in order to 
identify differentially expressed genes in the uterine 
endometrium on D12 and D15 of both the estrous cycle and 
of pregnancy. We compared microarray data from D12 and 
D15 of pregnancy to those from D12 and D15 of the estrous 
cycle and identified pregnancy-related variations in the 
patterns of global gene expression. 

MATERIALS AND METHODS 

Animals and tissue preparation 

All experimental procedures involving animals were 



conducted in accordance with the Guide for Care and Use 
of Research Animals in Teaching and Research and 
approved by the Institutional Animal Care and Use 
Committee of Yonsei University. Crossbred sexually 
matured gilts that were either not mated or were mated with 
fertile boars after the onset of estrus were assigned to the 
D12 or D15 estrous cycle groups or the D12 or D15 
pregnancy groups, respectively (n = 3 per d/status). 
Pregnancy was confirmed by the presence of an apparently 
normal conceptus in uterine flushings. Endometrium 
dissected from the myometrium was collected from the 
middle portion of the uterine horn. Endometrial tissues were 
snap-frozen in liquid nitrogen and stored at -80°C for RNA 
extraction. 

RNA isolation 

Total RNA from endometrial tissues was extracted using 
TRIzol (Invitrogen, Carlsbad, CA, USA) and the RNeasy 
Mini kit (Qiagen, Valencia, CA, USA) according to the 
manufacturer's instructions. The purity and integrity of the 
total RNA were checked using NanoDrop (NanoDrop 
Technologies, Wilmington, DE, USA) and Experion (Bio- 
rad, Hercules, CA, USA), respectively. 

Microarray analysis 

Five micrograms of total RNA were used for labeUng. 
Probe synthesis from the total RNA samples, hybridization, 
detection, and scanning were performed according to 
standard protocols from Affymetrix Inc. (Santa Clara, CA, 
USA). Briefly, cDNA was synthesized from total RNA 
using the One-Cycle cDNA Synthesis Kit (Affymetrix). 
Single-stranded cDNA was synthesized using Superscript II 
reverse transcriptase and T7-oUgo primers at 42°C for 1 h. 
Double-stranded (ds)-cDNA was obtained using DNA 
ligase, DNA polymerase I, and RNase H at 16°C for 2 h, 
followed by T4 DNA polymerase at 16°C for 5 min for gap 
filling. After clean-up with a Sample Cleanup Module 
(Affymetrix), ds-cDNA was used for in vitro transcription 
(IVT). cDNA was transcribed using the GeneChip IVT 
LabeUng Kit (Affymetrix) in the presence of biotin-labeled 
CTP and UTP. Ten micrograms of labeled cRNA were 
fragmented to 35-200 bp using fragmentation buffer 
(Affymetrix). Fragmented cRNA was hybridized to the 
porcine genome (Affymetrix) at 45 °C for 16 h according to 
the Affymetrix standard protocol. After hybridization, the 
arrays were washed in a GeneChip Fluidics Station 450 
with a non-stringent wash buffer at 25°C followed by a 
stringent wash buffer at 50°C. After washing, the arrays 
were stained with a streptavidin-phycoerythrin complex. 
After staining, intensities were determined with a GeneChip 
scanner 3000 (Affymetrix) controlled by GeneChip 
Operating Software (Affymetrix). 
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Data analysis 

Expression profiles were analyzed using GeneChip 
Operating Software (Affymetrix) to identify changes in 
expression levels in the uterine endometrium during the 
estrous cycle and early pregnancy. The operating software 
determined absolute analysis metrics (Detection, Detection 
p-value) using the scanned probe array data and compared 
the results between the different treatment group signals to 
generate the Change, Change p-value, and Signal log ratio 
(fold change) metrics after pre-processing of the raw data 
using the MASS algorithm. For normalization, data from 
each expression array were scaled to normalize the overall 
fluorescence intensity across each chip (average target 
intensity set at 500). The one-sided Wilcoxon's signed rank 
test was employed to generate the detection p-value (p< 
0.05). Two sets of algorithms were generated and used to 
determine change significance and change quantity metrics 
for every probe set. The change algorithm generated a 
Change p-value and an associated fold-change value, while 
the second algorithm gave a quantitative estimate of the 
change in gene expression in the form of the Signal log 
ratio. In the present study, the level of gene expression was 
considered to be increased if the Change p-value was less 
than 0.0025 and was considered to be decreased if the 
Change p-value was greater than 0.9975. Only relative 
changes equal to or greater than that of a two-fold level of 
expression were considered. We considered genes that 
matched the established Change p-value and the fold- 
change value to be differentially expressed genes. Probe 
identification was obtained using the gene ontology mining 
tool NetAffx (http://www.affymetrix.com/analysis/index.affx) 
and human homologues (Tsai et al., 2006). 

Functional annotation clustering analysis 

To evaluate biological functions of differentially 
expressed genes affected by pregnancy status in the uterine 
endometrium, we conducted functional annotation 
clustering analysis throughout DAVID program (Database 
for Annotation, Visualization, and Integrated Discovery, 
http://david.niaid.nih.gov/david/version2/index.htm) (Glynn 
et al., 2003). Official gene symbols of differentially 
expressed genes were submitted to DAVID for analysis of 
functional annotation clusters. Functional annotation 
clusters were determined by gathering related Gene 
Ontology (GO) terms with respect to biological process, 
cellular component, and molecular function. Because of 
limited number of annotated genes for porcine genome in 
NetAffx resulting from the Umited availabiUty of full-length 
porcine cDNA sequence, we used human gene symbols 
annotated for the Affymetrix porcine genome microarray 
probe identifications for the genes whose annotation was 
not available in NetAffx, as described in Tsai et al. (2006). 



Functional annotation clustering analysis using porcine and 
human gene symbols in DAVID program was conducted at 
highest stringency to clarify biological function of DEGs in 
the uterine endometrium. Functional groups were ordered 
by their overall enrichment score for the group based on the 
EASE scores of each term annotated. 

Quantitative real-time RT-PCR 

To analyze levels of selected differentially expressed 
gene mRNAs in the uterine endometrium on D12 and D15 
of the estrous cycle and pregnancy, real-time RT-PCR was 
performed using the Applied Biosystems StepOnePlus 
System (Applied Biosystems, Foster City, CA, USA) using 
the SYBR Green method. Complementary DNAs were 
synthesized from 4 \ig of total RNA isolated from different 
uterine endometrial tissues, and newly synthesized cDNAs 
(total volume of 21 |il) were diluted 1:4 with sterile water 
and then used for PCR. Specific primers based on porcine 
MUC5AC (GenBank accession number U10281.1; forward, 
5 '- CCA TCA TTT ACG AAG AGA CCG ACC-3 ' ; reverse, 
5'- GCC AGG TTT CAC CCT TCA TTC T-3'), NRIPl 
(GenBank accession number BV726910.1; forward, 5'- 
TCA AAA CTC CCC TGA GTC CTC CTT TC-3'; reverse, 
5'- ACG TCC TTG TCT TGT GTT TCT CGA CT-3'), 
RBP7 (GenBank accession number NM_00 1145222.1; 
forward, 5 '- GAG CAA AAT GGG GAT TCT TTT ACC A- 
3'; reverse, 5'- AC A CTT GGC CTT CAC AGA AC A TCT 
C-3'), ADAMTS20 (GenBank accession number 
XM_003 12662 1.1; forward, 5'- GTT GCT GAT GGT ACT 
CCT TGT GGA A-3'; reverse, 5'- CGA TAA CGC CAG 
GTA ATT GTC ATC T-3'), and porcine ribosomal protein 
L7 {RPL7), (GenBank accession number NM_001113217; 
forward, 5'-AAG CCA AGC ACT ATC ACA AGG AAT 
ACA-3'; reverse, 5'-TGC A AC ACC TTT CTG ACC TTT 
GG-3') were designed to amphfy cDNAs of 216 bp, 250 bp, 
250 bp, 285 bp, and 172 bp, respectively. The Power SYBR 
Green PCR Master Mix (Applied Biosystems) was used for 
PCR reactions. Final reaction volume was 20 |il including 2 
\i\ of cDNA, 10 ^1 of 2X Master mix, 2 ^1 of each primer, 
and 4 nl of dHzO. PCR conditions were 95°C for 15 min 
followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, and 
72°C for 30 s. The results are reported as the expression 
relative to the level detected on D12 of the estrous cycle 
after normalization of the transcript amount to the 
endogenous RPL7 control by the 2"'^''^ method (Livak and 
Schmittgen, 2001). 

Statistic analysis 

Data from real-time RT-PCR analysis for comparison of 
expression levels of selected genes in endometrium on D12 
and D15 of the estrous cycle and pregnancy were subjected 
to the Student's t test procedure of SAS, and are presented 
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as means with standard errors. 

RESULTS 

Comparison of gene expression in the uterine 
endometrium between D12 and D15 of pregnancy and 
D12 and D15 of the estrous cycle 

To determine global gene expression profiles in the 
uterine endometrium on D12 and D15 of the estrous cycle 
and of pregnancy in pigs, microarray analyses were 
performed using a total of 12 Affymetrix porcine genome 
microarray chips. To understand pregnancy-specific 
changes in endometrial gene expression, we identified 
differentially expressed genes with at least a two -fold 
increase or decrease in the uterine endometrium between 
D12 of pregnancy and D12 of the estrous cycle, as well as 
those with at least a two-fold increase or decrease between 
D15 of pregnancy and D15 of the estrous cycle. 

After analyzing the raw data using the MAS5 algorithm 
of the GeneChip Operating Software, differentially 
expressed genes among the compared pairs were obtained. 
We considered genes that were below a Detection p-value 
of 0.05 and a two-fold change to be differentially expressed 
genes. Among 23,937 genes, 99 and 213 genes were up- 
regulated in the uterine endometrium on D12 and D15 of 
pregnancy compared to D12 and D15 of the estrous cycle, 
respectively. In addition, 92 and 231 genes were down- 
regulated in the endometrium on D12 and D15 of 
pregnancy compared to D12 and D15 of the estrous cycle, 
respectively. Volcano plot images showing comparative 
analyses of endometrial gene expression between the 
estrous cycle and pregnancy are shown in Figure 1 . 

Genes differentially expressed on D12 of pregnancy 



compared to those of D12 of the estrous cycle 

Genes that had higher expression in the uterine 
endometrium on D12 of pregnancy compared to those of 
D12 of the estrous cycle included sulfotransferase family 
IE (SULTlEl), LPAR3, caspase 3 {CASP3), interleukin 1 
receptor accessory protein (ILIRAP), and fibroblast growth 
factor 7 (FGFT) (Supplementary Table 1). Genes that had 
lower expression in the endometrium on D12 of pregnancy 
compared to those of D12 of the estrous cycle included 
glutathione S-transferase (GST), mucin (MUC5AC), solute 
can-ier family 2 {SLC2A2), cytochrome P450 3A46 
(CYP3A46), keratin 7 (KRT7), chemokine ligand 12 
(CXCL12), and interferon regulatory factor 7 (IRF7) 
(Supplementary Table 2). 

Genes differentially expressed on D15 of pregnancy 
compared to those of D15 of the estrous cycle 

Genes that had higher expression in the uterine 
endometrium on D15 of pregnancy compared to those of 
D15 of the estrous cycle included endoplasmic reticulum 
degradation enhancer mannosidase alpha-like 3 {EDEM3), 
neurotrimin (NTM), UDP-glucose ceramide glucosyltransferase 
(UCGC), nuclear receptor interacting protein 1 (NRIPl), 
and SMAD family member 2 (SMAD2) (Supplementary 
Table 3). Genes that had lower expression in the 
endometrium on D15 of pregnancy compared to those of 
D15 of the estrous cycle included cellular retinol binding 
protein 7 {RBP7), parathyroid hormone (PTH), and myosin 
VI (MY06) (Supplementary Table 4). 

Functional annotation clustering analysis 

To identify putative biological functions regulated by 
pregnancy in the uterine endometrium, we performed 
functional annotation clustering analysis of differentially 




Figure 1. Volcano plot images for the comparative analysis of endometrial gene expression on day (D) 12 of pregnancy to D12 of the 
estrous cycle (A) and on D15 of pregnancy to D15 of the estrous cycle (B) in pigs. The X axis represents the average log_2 (fold- 
change), and the Y axis represents the -log (p-value) with base 10. Location in the upper left shaded area implies down-regulation; 
location in the upper right shaded area implies up-regulation. Cutoff values were below p<0.05 (horizontal line) and above an average 
two-fold change (vertical lines). 
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Figure 2. Functional annotation ciiarts siiowing distribution of ontological analysis of up-regulated (A) and down-regulated genes (B) in 
the uterine endometrium on day (D) 12 of pregnancy compared with those on D12 of the estrous cycle, and up-regulated (C) and down- 
regulated genes (D) on D15 of pregnancy compared with those on D15 of the estrous cycle. Biological terms were identified by genes 
with respect to biological process at level one through DAVID. Gene ontology terms or biological features of differently expressed genes 
and the number of genes in each category are indicated. 



expressed genes by the DAVID program. As a result, we 
obtained enriched functional groups in each cluster and 
clarified gene ontology (GO) terms or biological features of 
DEGs in the uterine endometrium (Figure 2). In genes up- 
regulated in the uterine endometrium on D12 of pregnancy 
compared with those on D 12 of the estrous cycle, functional 
annotation clustering analysis showed that functional 
groups of membrane fraction and endoplasmic reticulum 
were enriched biological feature (Supplementary Table 5). 
Biological functions of cell adhesion and immune cell 
activation were enriched in genes down-regulated in uterine 
endometrium on D 12 of pregnancy compared with those on 
D12 of the estrous cycle (Supplementary Table 6). RNA 
splicing and mRNA metabolism were enriched in genes up- 
regulated in the uterine endometrium on D15 of pregnancy 
compared with those on D15 of the estrous cycle 
(Supplementary Table 7), whereas cytoskeloton, micro- 



tubule, actin dynein related functions were enriched in 
genes down-regulated in the uterine endometrium on D15 
of pregnancy compared to those on D15 of the estrous cycle 
(Supplementary Table 8). 

Real-time RT-PCR analysis of ADAMTS20, MUC5AC, 
NRIPl, and RBP7 in the uterine endometria on D12 and 
D15 of the estrous cycle and pregnancy 

To determine validity of our microanalysis data, we 
selected ADAMTS20, MUC5AC, NRIPl, and RBP7 from 
genes up-regulated or down-regulated in comparison groups 
between D12 and D15 of the estrous cycle and pregnancy, 
and analyzed expression levels of those genes using real- 
time RT-PCR analysis. Levels of ADAMTS20 expression 
were significantly higher in the uterine endometrium on 
D12 of pregnancy compared with those on D12 of the 
estrous cycle (p<0.05). Expression levels of MUC5AC were 
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significantly lower in the uterine endometrium on D12 of 
pregnancy compared with those on D12 of the estrous cycle 
(p<0.05). Levels of NRIPl expression was significantly 
higher in the uterine endometrium on D15 of the pregnancy 
compared with those on D15 of the estrous cycle (p<0.05). 
However, levels of RBP7 expression were not different in 
the uterine endometrium between D15 of the estrous cycle 
and pregnancy (p>0.05). 

DISCUSSION 

The estabUshment of pregnancy is a complex process 
that requires well-coordinated interactions between the 
maternal uterus and the developing conceptus. In pigs, the 
implantation process begins around D12 of pregnancy and 
completes around D18 (Geisert and YeUch, 1997). During 
this period, the presence of the conceptus in the uterus 
induces changes in endometrial gene expression in 
preparation for conceptus attachment to the endometrial 
epithelial cells and extension of the life span of the corpus 
luteum in the ovaries (Spencer and Bazer, 2004). To better 
understand the implantation process, this study aimed to 
identify changes in endometrial gene expression on D12 
and D15 of pregnancy in comparison to expression on D12 
and D15 of the estrous cycle using a microarray technique. 
In this study, we identified genes that were up-regulated or 
down -regulated in the uterine endometrium during the 
implantation period. These genes that were differentially 
expressed according to pregnancy status included genes 
involved in immunity, steroidogenesis, cell-to-cell 
interaction, and tissue remodeling. Furthermore, functional 
annotation functional clustering analysis of these genes 
showed that biological functions of ceU-to-cell interaction 
and immune cell activation were decreased in the uterine 
endometrium on D12 of pregnancy compared to those on 
D12 of the estrous cycle and RN A processing was increased 
in the uterine endometrium on D 15 of pregnancy compared 
with those on D15 of the estrous cycle. 

On around D12 of pregnancy in pigs, maternal 
recognition of pregnancy occurs through a mechanism 
explained by the endocrine-exocrine theory (Bazer and 
Thatcher, 1977). The elongating conceptus secretes estrogen 
into the uterine lumen, which, in turn, acts on the uterine 
endometrium and redirects the secretion of PGF2a toward 
the uterine lumen (exocrine) so that the corpus luteum (CL) 
function is maintained in the ovary. On the other hand, the 
absence of a conceptus in cyclic pigs on D12 causes 
secretion of PGF2a into the uterine vasculature (endocrine), 
which exerts luteolytic action on the CL. Thus, it is likely 
that the expression of endometrial genes on D12 of 
pregnancy is affected by the action of the estrogen 
originating from the conceptus. Indeed, it has been shown 
that the expressions of FGF7 and LPAR3 are up-regulated 



following estrogen treatment in uterine endometrium (Ka et 
al., 2008; Seo et al., 2008). In addition to estrogen, the 
conceptus also produces cytokines such as IL-ip, IFN-y and 
IFN-5 (Spencer and Bazer, 2004), and the uterine 
environment during this period is also influenced by 
progesterone of ovarian origin. Therefore, expression of 
endometrial genes on D12 of pregnancy would be affected 
by a coordinated interaction of these factors. 

Genes up-regulated in the endometrium of pregnant pigs 
on D12 of pregnancy compared to genes of D12 of the 
estrous cycle include SULTlEl, LPAR3, ILIRAP, and 
FGF7. Previous studies in pigs have demonstrated that 
these genes are expressed in the endometrium at 
significantly higher levels on D12 of pregnancy than on 
D12 of the estrous cycle (Ka et al., 2000; Kim et al., 2002; 
Ross et al., 2003b; Seo et al., 2008), indicating the validity 
of our microarray approach in analyzing gene expression. In 
addition, many other genes, including ADAMTS17, 
ADAMTS20, CASP3, CCNA2, and CYP19, were also 
identified as being up-regulated in the endometrium of 
pregnant pigs on D 12 of pregnancy compared with genes on 
D12 of the estrous cycle. In functional annotation clustering 
analysis using genes up-regulated in the uterine 
endometrium on D 12 of pregnancy compared with genes on 
D12 of the estrous cycle, any notable biological function 
was not found due to low relatedness of biological terms 
and low numbers of genes having similar function. 

Down-regulated genes in the endometrium of pigs on 
D12 of pregnancy compared with genes on D12 of the 
estrous cycle include genes representing biological 
functions of leukocyte activation (CDS6, CD8A, CXCL12, 
KLRKl, MLL5, PTPN22, and UNC13D), cell adhesion 
(AJAPl, CNTN5, CXCL12, CYFIP2, EGFL6, FATl, ITGB4, 
ITGB7, MUC5AC, PKP2, and ROBOl), and cell junction 
(AJAPl, CYFIP2, GABRA5, GRIPl, PKP2, SHR00M2, 
and TNS3). During the implantation period in pigs, immune 
cell population changes and effector T cell population 
decreases in the subepithelial regions of the uterine 
endometrium (Croy et al., 1987; Bischof et al., 1996; 
Kaeoket et al., 2003), suggesting that maternal immune 
suppression contributes to the establishment of pregnancy. 
CD86, together with CD80, regulates T cell function as a 
ligand for CD28 and cytotoxic lymphocyte antigen-4 
(CTLA-4), and expression of CD86 is detected in 
macrophages and dendritic cells at the feto-maternal 
interface in human (Banchereau et al., 2000) and mouse 
(Repnik et al., 2008). Interestingly, decreased expression of 
CD86 in the uterine endometrium with the conceptus on 
D12 of pregnancy compared D12 of the estrous cycle 
indicate that regulation of maternal uterine T cell function 
may also be important for the establishment of pregnancy in 
pigs. CXCL12 plays a key role in lymphocyte trafficking, 
cellular proUferation, organogenesis, vascularization, and 
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Figure 3. Real-time RT-PCR analysis of selected differentially expressed genes, ADAMTS20 (A), MUC5AC (B), NRIPl (C), and RBP7 
(D), in the uterine endometrium during the implantation period in pigs. Levels of ADAMTS20 and MUC5AC mRNA were significantly 
different in the endometrium on D12 of the estrous cycle and pregnancy (p<0.05). Levels of NRIPl mRNA were significantly higher in 
the endometrium on D15 of pregnancy than D15 of the estrous cycle (p<0.05), but levels of RBP7 were not different (p>0.05). 
Abundance of mRNA is presented as the expression relative to the levels of ADAMTS20, MUC5AC, NRIPl and RBP7 mRNAs measured 
in the uterine endometrium after normaUzation of the transcript amount to RPL7 mRNA. Data are presented as means with standard 
error. 



embryogenesis (Murdoch, 2000; Wu et al., 2004). CXCL12 
has been shown to be expressed in the uterine endometrium 
of humans and bovines (Hanna et al., 2003; Mansouri-Attia 
et al., 2009). In human trophoblasts and decidual cells, 
CXCL12 increases the invasiveness of trophoblasts and 
matrix metalloproteinase (MMP) 9 and MMP2 activity 
(Zhou et al., 2008) and has also been shown to be 
associated with endovascular invasion of CD16-NK cells in 
hviman placenta (Hanna et al., 2003). In this study, 
expression of CXCL12 decreased in pregnant uterine 
endometrium. Although its function in the uterine 
endometrium needs to be further elucidated during early 
pregnancy in pigs, modulation of CXCL12 action may also 
be important for embryo implantation. 

During the implantation period in rodents, humans, and 
ruminants, it has been shown that expression of MUCl, a 
type of transmembrane mucin, is decreased in uterine 
endometrial epithelial cells to allow for adhesion of the 
conceptus trophectoderm to endometrial epithelial cells 
(Carson et al., 1998; Spencer et al., 2004). MUC5AC is a 
secreted, large, gel-forming type of mucin (Gendler and 
Spicer, 1995). Its expression has been measured in human 
uterine endocervical epithelial cells as a major form of gel- 
forming mucin in the menstrual cycle (Gipson et al., 1997), 



and it is associated with a change in eosinophilic cells in 
endometrial carcinomas (Moritani et al., 2005). Decreased 
expression of MUC5AC in pregnant uterine endometrium 
suggests that MUC5AC, in addition to MUCl, may block 
adhesion of the conceptus trophectoderm and must be 
removed for the trophectoderm to undergo epithelial cell 
attachment. 

We also compared the genes expressed in the 
endometrium on D15 of pregnancy to those expressed on 
D15 of the estrous cycle. On D15 of the estrous cycle, the 
period of late diestrus in pigs, luteolysis is initiated by 
extended exposure of the endometrium to progesterone for 
10 to 12 d, followed by increased pulsatile production of 
PGF2a (Ziecik, 2002). According to a study in sheep, the 
initial exposure to progesterone down-regulates the 
endometrial epithelial expressions of ESRl and PGR, and 
the loss of the action of progesterone permits estrogen to 
increase ESRl and OXTR in the uterine endometrium 
(Spencer and Bazer, 2004). Pulsatile release of oxytocin 
induces subluteolytic pulses of PGFaa from the uterine 
endometrium. PGF2a causes a supplemental release of 
oxytocin from the CL, which in turn stimulates the 
secretion of a luteolytic pulse of PGF2a from the uterine 
endometrium through a positive feedback mechanism. This 
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positive feedback loop continues until the CL is depleted of 
oxytocin (Spencer and Bazer, 2004). On the other hand, on 
D15 of pregnancy, maternal recognition of pregnancy is 
established, and interdigitation between the chorion and 
endometrial epithelial cells develops to form the 
epitheliochorial placenta (Dantzer, 1985). During this time, 
the uterine environment is under the influence of 
progesterone, which is continuously produced from the 
ovaries, and cytokines originating from the conceptus, such 
as IFN-y and IFN-5 (Cencic et al., 2003). It has been shown 
that a large amount of IFN-y is secreted by the pig 
conceptus between d 12 and 20 of pregnancy (up to 250 mg 
per uterine horn), with the highest levels of synthesis on 
D15 and DI6 (La Bonnardiere et al., 1991). At this time, 
ESRl expression is decreased in the endometrium, and 
PGR is only detectable in stromal cells (Geisert et al., 1993; 
1994). 

In this study, many novel genes whose roles in the 
uterus have not been fully studied were identified as being 
differentially expressed in the endometrium between D15 of 
pregnancy and D15 of the estrous cycle. Among genes that 
were up-regulated in the endometrium on D15 of pregnancy 
compared to those on D15 of the estrous cycle were AKAPl, 
AKAPll, AKAP13, CLCN3, EDEM3, NRIPl, and SMAD2; 
down-regulated genes in the endometrium on D15 of 
pregnancy compared to those on D15 of the estrous cycle 
included ATP2C2, EFHA2, F0XN3, MY06, PERI, PTH 
and RBP7. NRIPl, formerly known as receptor interacting 
protein 140 (RIP140), is a widely expressed corepressor 
that has the potential to inhibit the transcriptional activities 
of most nuclear receptors, including ESRl, PGR, and GR. 
Studies have shown that NRIPl is expressed in metabolic 
tissues and plays a critical role in lipid metabolism (Fritah, 
2009). The importance of NRIPl in reproduction has been 
shown in a study of Nripl-mxM mice (White et al., 2000), 
which are viable, but whose females are infertile due to a 
failure to release mature oocytes (White et al., 2000). The 
expression and function of NRIPl in reproductive tissues 
have not been heavily studied, but some data in mice and 
rats have shown that Nripl is expressed in the uterus 
(Nephew et al., 2000; Leonardsson et al., 2002). In mice, 
Nripl expression is localized to glandular epithelial cells, 
stromal cells, and decidual cells, but is not observed in 
luminal epithelial cells in the uterine endometrium 
(Leonardsson et al., 2002). In women, mutation of the 
NRIPl gene is associated with endometriosis, suggesting 
that mutation of this gene may be a predisposing factor for 
endometriosis (Caballero et al., 2005). There is no 
information available on the expression and function of 
NRIPl in the porcine uterus. However, NRIPl may play an 
important role in regulating lipid hormonal actions for the 
establishment of pregnancy, because the uterus is 
continuously influenced by steroid hormones (mainly 



estrogen and progesterone) during early pregnancy. Further 
characterization of the mechanism of NRIPl action in the 
uterus during early pregnancy is needed. 

SMAD2 is one of the key modulators in the TGF-P 
signaling pathway. When the TGF receptor is activated by 
TGF-P binding, SMAD2 binds to the TGF-P-receptor 
complex and is activated by C-terminal phosphorylation. 
Phosphorylated SMAD2 pairs with SMAD4 and relocates 
to the nucleus, where it acts as a transcription factor to 
induce target gene expression (ten Dijke and Hill, 2004). In 
rodents, SMAD2 expression in the uterus has been well 
characterized. During the estrous cycle, Smad2 expression 
is locahzed to luminal and glandular epitheUal cells, but its 
expression is predominantly localized to the subluminal 
stroma surrounding the implanting embryo during the peri- 
implantation period (Liu et al., 2004). In the porcine uterus, 
phosphorylated SMAD2/3 was localized to endometrial 
luminal and glandular epithelilal cells, stromal cells, and 
endothelial cells (as well as the conceptus trophectoderm) 
on D13 of pregnancy, and activation of SMAD2/3 in 
endometrial luminal epithelial cells by TGF-P infusion into 
the uterine lumen has been demonstrated (Massuto et al., 
2010). In pigs, it has been suggested that TGF-p, which is 
present at the maternal-fetal interface during the 
implantation period (Gupta et al., 1998), acts in 
trophectodermal cell-to-luminal epithelial cell adhesion 
(Massuto et al., 2010). Thus, increased expression of 
SMAD2 in the endometrium of pregnant pigs on D15 may 
indicate that TGF-P signaUng is activated in conceptus 
implantation. 

Recently, microarray analysis has been appUed in 
determining DEGs in the uterine endometrium in pregnant 
pigs compared with genes expressed in non-pregnant pigs, 
showing that DEGs are involved in gene function of 
developmental process, transporter activity, calcium ion 
binding, apoptosis, cell motility, enzyme-linked protein 
receptor signaling pathway, positive regulation of cell 
proUferation, ion homeostasis, and hormone activity 
(0strup et al., 2010). Microarray analysis has also been 
applied to determine the pattern of global gene expression 
during the implantation period in several species, including 
humans (Kao et al., 2002), mice (Reese et al., 2001), bovine 
(Ushizawa et al., 2004), and ovine (Satterfield et al., 2009). 
A direct comparison of genes expressed in the uterine 
endometrium during the implantation period in pigs to 
genes in other species is not possible because of varying 
experimental settings between studies and hmited 
information on annotated genes in the porcine genome. 
However, it seems that very few genes expressed in the 
porcine uterine endometrium match the genes expressed in 
other species, which may result from different mechanisms 
of implantation and placentation in pigs. 
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The data presented in this study provides information on 
the genes expressed in the uterine endometrium during the 
implantation period in pigs. Expression of a gene during 
this period does not necessarily mean that the gene is 
critical for conceptus implantation; as such, further 
functional characterization of those genes is essential. 
Nevertheless, data from this study will provide insight into 
the mechanism of the implantation process for the 
establishment of pregnancy in pigs. In addition, the 
experimental setting used in this study may be applicable to 
the screening of candidate genes related to prolificacy 
among pig breeds and to low implantation rates in the 
production of somatic cell nuclear transfer cloned pigs. 
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in the uterine endometrium on day (D) 12 of pregnancy compared to ^ 
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Probe sel ID 


.'\verage 
k)g_2(r()ld change; 


Ciene sMiibol 


Cleiie lille 


RefSeq 
tiaiiscripl ID 


Ssc.l2.1.Sl_at 


5.1138 


SULTlEl 


Sulfotransferase family IE, estrogen-preferring, member 1 


NM_213992 


Ssc.31161.1.Al_at 


4.2407 


ADAMTS20 
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Ubiquitin protein hgase E3 component n-recognin 7 (putative) 


NM_001 100417 


Ssc.l2523.1.Al_at 


2.3272 


CLIC4 


Chloride intracellular channel 4 


M"\ 4 m '2r\A'j 
JNJVl_01 jy4j 


Ssc.l5886.1.Sl_at 


2.3001 


CASP3 


Caspase 3, apoptosis-related cysteine peptidase 


M\ 4 O 1 /I 1 1 1 


Ssc.7864.1.Al_at 


2.1904 


ILIRAP 


Interleukin 1 receptor accessory protein 


iMVl_UUZ i oZ 


Ssc.l8027.1.SI_at 


2.1855 


BTBD3 


BTB (POZ) domain containing 3 


NM_0 14962 


Ssc.6441.2.Sl_at 


2.1576 


DHRS7 


Dehydrogenase/reductase (SDR family) member 7 


NM_0 16029 


Ssc.8552.3.Sl_a_at 


2.1220 


GGTAl 


Glycoprotein, alpha-galactosyltransferase 1 


NM_213810 


Ssc.21460.I.AI_at 


1 .9430 


CREGl 


Cellular repressor of ElA-stimulated genes 1 


NM_003851 


Ssc.l3221.I.Al_at 


1.8729 


PC 


Pyruvate carboxylase 


NM_000920 


Ssc.4999.1.SI_at 


1.8630 


AVPII 


Arginine vasopressin-induced 1 


NM_021732 


Ssc.l4304.1.Al_at 


1.8514 


CLSTN2 


Calsyntenin 2 


NM_022131 


Ssc.l5923.1.Al_at 


1.8044 


FGF7 


Fibroblast growth factor 7 


AF2 17463.1 


Ssc.30862.1.SI_at 


1.8001 


DNAJB9 


DnaJ (Hsp40) homolog, subfamily B, member 9 


NM_0 12328 


Ssc.31013.2.SI_at 


1.7794 


KCTD3 


Potassium channel tetramerisation domain containing 3 


NM_016121 


Ssc. 12056. 1.AI_at 


1.7746 


L0XL4 


Lysyl oxidase-like 4 


NM_032211 


Ssc.24450.I.Sl_at 


1.7716 


F13B 


Coagulation factor XIII, B polypeptide 


NM_00I994 


Ssc.l2412.1.Al_at 


1.7315 


SDF2L1 


Stromal cell-derived factor 2-like 1 


NM_022044 


Table 2. List of top 30 ; 


genes down-regulated in the uterine endometrium on day (D) 12 of pregnancy compared to i 


genes on D12 of the 


estrous cycle by average fold change 








Probe set ID 


Average 
log_2(Fold change) 


Gene symbol 


Gene title 


RefSeq 
transciipt ID 


Ssc.l0837.1.Al_at 


-3.994277 


ROBOl 


Roimdabout, axon guidance receptor, homolog 1 (Drosophila) 


NM_002941 


Ssc.I5850.1.Sl_a_at 


-3.892519 


TRAm 


T cell receptor alpha locus 


NW_001838110.1 


Ssc.I6377.2.Al_at 


-3.697358 


GSTA3 


Glutathione S-transferase alpha 3 


NM_000847 


Ssc.3904.1.SI_at 


-2.781738 


RPS6KA6 


Ribosomal protein S6 kinase, 90kDa, polypeptide 6 


NM_014496 


Ssc.l6767.LSl_at 


-2.76431 


HNRNPL 


Heterogeneous nuclear ribonucleoprotein L-like 


NM_001005335 


Ssc.338.1.Sl_at 


-2.656394 


MUC5AC 


Mucin 5AC, oligomeric mucus/gel-forming 


NM_017511.1 


Ssc.22190.I.Sl_at 


-2.632552 


AJAPl 


Adherens junctions associated protein 1 


NM_001042478 


Ssc.8621.1.Al_at 


-2.490124 


PRKACB 


Protein kinase. cAMP-dcpcndent. catalytic. Beta 


NM_002731 


Ssc.23849.1.AI_at 


-2.457561 


SLC2A2 


Solute carrier family 2 (facilitated glucose transporter), member 2 


NM_000340 


Ssc.25254.1.Sl_at 


-2.378661 


RHOJ 


Ras homolog gene family, member I 


JNM_0z06OJ 


Ssc.21224.1.Al_at 


-2.352619 


KRT7 


Keratin 7 


NM_005556 


Ssc.6556.2.Al_at 


-2.344349 


GABRA5 


Gamma-aminobutyric acid (GABA) A receptor, alpha 5 


NM_000810 


Ssc.27912.1.Sl_at 


-2.205287 


ZNF608 


Zinc finger protein 608 


NM_020747 


Ssc.24434.1.Sl_at 


-2.16435 


C14orfl45 


Chromosome 14 open reading frame 145 


NM_152446.3 


Ssc.26326.1.Al_at 


-2.134067 


CYP3A46 


Cytochrome P450 3A46 


NM_001134824 


Ssc.ll260.1.Al_at 


-2.124412 


GSTAl 


Glutathione S-transferase alpha 1 


NM_145740 


Ssc.8053.1.Sl_at 


-1.953481 


POLL 


Polymerase (DNA directed), lambda 


NM_013274 


Ssc.2618.1.Sl_at 


-1.827725 


MAMDC4 


MAM domain containing 4 


NM_206920 


Ssc.5641.1.Sl_at 


-1.80541 


ANK2 


Ankyrin 2, netironal 


NM_001127493 


Ssc.l8992.1.Al_at 


-1.784956 


ANKSD6 


Ankyrin repeat domain 6 


NM_014942 


Ssc.30445.1.Al_at 


-1.761477 


ALK 


Anaplastic lymphoma receptor tyrosine Kinase 


NM_004304 


Ssc.27318.1.Sl_at 


-1.757688 


KRT7 


Keratin 7 


NM_005556 


Ssc.30532.1.Al_at 


-1.725776 


XRCC2 


X-ray repair complementing defective repair in Chinese hamster 

cells 2 

EGF-like-domain. multiple 6 


NM_005431 


Ssc.27626.I.Sl_at 


-1.719731 


EGFL6 


NM_015507 


Ssc. 12223. 1.AI_at 


-1.690582 


PTPN22 


Protein tyrosine phosphatase, non-receptor type 22 (lymphoid) 


NM_012411 


Ssc.336.1.Sl_at 


-1.631299 


USPIS 


Ubiquitin specific peptidase 1 8 


NM_2 13826 


Ssc.l5708.1.Sl_at 


-1.621459 


HBZ 


Hemoglobin, zeta 


NMJX)5332 


Ssc.7243.1.Al_at 


-1.619644 


CXCLI2 


Chemokine (C-X-C motif) ligand 12 (stromal cell-derived factor 1) 


NMJX) 1009580 


Ssc.21162.1.Sl_s_at 


-1.609334 


IRF7 


Interferon regulatory factor 7 


NMJX) 1097428 


Ssc.286.1.Sl_s_at 


-1.553267 


IRG6 


Inflarmnatory response protein 6 


NM_213817 
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Table 3. List of top 30 


genes up-regulated in the uterine endometrium on day (D) 15 of pregnancy compared to \ 


jenes on D15 of the 


estrous cycle by average fold change 










Average 
log_2(Fold change) 


Gene symbol 


Gene title 


RefSeq 
transcript ID 


Ssc.28494.1.Sl_at 


3.9383217 


RAB5B 


RAB5B, member RAS oncogene family 




Ssc.I8273.1.Al_at 


3.9374428 


MAGIl 


Membrane associated guanylate kinase, WW and PDZ domain 
containing 1 


MM 17t'%l'% 

I1iVl_l 1 JJIJ 


Ssc.22734.1.Al_at 


3.505833 


EDEM3 


ER degradation enhancer, mannosidase alpha- like 3 


NM 025191 


Ssc.6392.1.Sl_at 


3.3603669 


NTM 


Neurotrimin 


NM 001144058.1 


Ssc.25740.1.Al_at 


3.3474057 


AKAPl 


A kinase (PRKA) anchor protein 1 


NM_003488 


Ssc.9329.1.Al_at 


3.0693123 


PTPLAD2 


Protein tyrosine phosphatase-like A domain containing 2 


MA/I nn 1 A 1 HQ 1 ^ 

IN ivi_uu 1 u 1 uy 1 J 


Ssc.l0593.2.Al_at 


3.0355328 


IFI44L 


Interferon-induced protein 44-like 


IN 1V1_1)V.'0 OZ.\) 


Ssc.ll925.I.AI_at 


2.8589389 


SELIL 


Sel-1 suppressor of lin-12-like (C. Elegans) 


IN1V1_UUjUDj 


Ssc.25678.I.Sl_at 


2.8213114 


SYTl 


Synaptotagmin I 


IN IV1_UU J vjy 


Ssc.l2809.5.Al_at 


2.7832047 


HBA2 


Hemoglobin, alpha 2; hemoglobin, alpha 1 


INIV1_UUUD1 / 


Ssc.7455.1.Al_at 


2.7595911 


AASDHPPT 


Aminoadipatc-semialdehydc dehydrogenase- 
phosphopantcthcinyl transferase 


IN iVl_U Lj'fZj 


Ssc.7682.1.AI_at 


2.6842858 


0SBPL9 


Oxystcrol binding protein-like 9 


IN oD 


AFFX-Ss_18SrRNA_at 


2.6236143 


AFFX-Ss_ 
ISSrRNA 


AFFX-Ss_ ISSrRNA 


Null 


Ssc.21585.2.Sl_at 


2.5961034 


UBR7 


Ubiquitin protein ligase E3 component n-reeognin 7 (putative) 


NM no 1 1 nod. 1 7 


Ssc.l3473.1.Al_at 


2.5538991 


UGCG 


UDP-glueose eeramide ghicosyltransferase 


NM 003358 


Ssc.27407.1.Al_at 


2.5424531 


NRIPI 


Nuclear receptor interacting protein 1 


N!M no'^dSQ 

INIVl \J\JJHC)y 


Ssc.3975.2.Al_at 


2.4775446 


STK38L 


Serine/threonine kinase 38 like 




Ssc.28192.1.Al_at 


2.4555041 


FAM153A 


Family with sequence similarity 153, member A 




Ssc.20353.1.Sl_at 


2.3855645 


MALTl 


Metastasis associated lung adenocarcinoma transcript 1 (non- 
protein coding) 


NR 002819.2 


Ssc.ll870.1.Al_at 


2.2433352 




Nuclear receptor interacting protein 1 


NM 003489 


Ssc.l7508.2.Al_at 


2.2266604 


RCAN2 


Regulator of calcineurin 2 


NM_005822 


Cop 814^ 0 Al flt 


2.1943771 


RAPGEF2 


Rap guanine nucleotide exchange factor (GEF) 2 


NM_014247 


Ssc.l7152.2.Al_at 


2.1849978 


C0Q2 


Coenzyme Q2 homolog, prenyltransferase (Yeast) 


NM_0I5697 


Ssc.27130.1.Al_at 


2.1615496 


SLC25A36 


Solute carrier family 25, member 36 


NM_001 104647 


Ssc.3325.3.Sl_at 


2.1184436 


EXTl 


Exostoses (multiple) 1 


NM_000127 


Ssc.22645.1.Sl_at 


2.0685056 


DDX46 


DEAD (Asp-Glu-Ala-Asp) box polypeptide 46 


NM 014829 


Ssc.6175.1.Sl_at 


2.0623846 


LCOR 


Ligand dependent nuclear receptor corepressor 


NM_032440 


Ssc.l3601.1.Al_at 


1.8954884 


MYOIB 


Myosin IB 


NM_001130158 


Ssc.24402.1.Sl_at 


1.8920959 


CCDC76 


Coiled-coil domain containmg 76 


NM_0 19083 


Ssc.27133.1.Al_at 


1.8789606 


TMEM68 


Transmembrane protein 68 


NM_152417 


Table 4. List of top 30 j 


jenes down-regulated in the uterine endometrium on day (D) 15 of pregnancy compared to 


genes on D15 of the 


estrous cycle by average fold change 








Probe sel ID 


Average 
log_2(rc)ld changej 


Ciene symbol 


(ieiie lille 


RefSeq 
liansciipl ID 


Ssc.25087.1.Sl_at 


-4.665941 


DNAJC13 


DnaJ (Hsp40) homolog, subfamily C, member 13 


NM 015268 


Ssc.2015.2.Sl_at 


^.469763 


ZNF449 


Zinc finger protein 449 


NM 152695 


Ssc.7817.1.Al_at 


^.302148 


NETOl 


Neuropilm (NRP) and toUoid (TIX)-like 1 


NM 138966 


Ssc.23424.1.Al_at 


^.072875 


DTNA 


Dystrobrevin, alpha 


NM_001128175 


Ssc.8989.1.Al_at 


-3.882042 


PPHLNl 


Periphihn 1 


NM 016488 


Ssc.21088.2.Sl_at 


-3.009555 


CSTFl 


Cleavage stimulation factor, 3' pre-RNA, subunit 1, 50kDa 


IN iyi_\J\Ji\JjDJi^i 


Ssc.9834.2.Sl_at 


-2.875995 


Null 


Sn36_F05.f sn Sus scrofa cDNA 


^^DH / 1 jXJjL. 1 


Ssc.l600.1.Al_at 


-2.684489 


RIMKLB 


Ribosomal modification protein rimK-like family member B 


IN IYl_U.iU 1 D't.Z, 


Ssc.l600.1.Al_a_at 


-2.616954 


RIMKLB 


Ribosomal modification protein rimK-like family member B 


NM 020734.2 


Ssc.27385.1.Sl_at 


-2.580643 


RBP7 


Retinol binding protein 7, cellular 


INIVl UUl iHJ^Zf^ 


Ssc.668.1.Sl_at 


-2.543521 


PTH 


Parathyroid hormone 


NM 91 Aim 

IN IVl 1 H-H-U 1 


Ssc.958.1.Sl_at 


-2.543376 


ALDOB 


Aldolase B, fruetose-bisphosphatc 


NM nono'^s 


Ssc.27953.1.Sl_at 


-2.289442 


MLL 


Myeloid/lymphoid or mixed-lineage leukemia (trithorax 
homolog, Drosophila) 


IN IVl I J\7 J 7J3 


Ssc.21468.2.Al^at 


-2.194019 


LARPl 


La ribonueleoprotein domain family, member 1 


NM 015315 


Ssc.30456.1.Al_at 


-2.155805 


CCDC135 


Coiled-coil domain containing 135 


NM 032269 


Ssc.27345.1.Sl_at 


-2.107692 


Null 


940760 MARC 4PIG Sus scrofa cDNA 


CN153625.1 


Ssc.7219.1.AI_at 


-1.926105 


ERC2 


ELKS/RAB6-interacting/CAST family member 2 


NM 015576 


Ssc.2279I.IAI_at 


-1.922899 


EE A I 


Early endosome antigen 1 


NM 003566 


Ssc.7672.I.AI_at 


-1.847485 


CAPZB 


Capping protein (actin filament) muscle Z-line, beta 


NM 004930 


Ssc.25130.1.Al_at 


-1.83943 


C6orf9S 


Chromosome 6 open reading frame 98 


NM 001099267 


Ssc.30379.1.Al_at 


-1.831398 


ARMC3 


Armadillo repeat containing 3 


NM_ 173081 


Ssc.7191.1.Al_at 


-1.805 


LSM14A 


LSM14A, SCD6 homolog A (S. Cerevisiae) 


NM_001 114093 


Ssc.2589.1.Sl_at 


-1.790171 


TNNI3 


Troponin 1 type 3 (cardiac) 


NM_001098599 


Ssc.23948.1.Al_at 


-1.727051 


DYNC1I2 


Dynein, cytoplasmic 1, intermediate chain 2 


NM_001378 


Ssc.8133.1.Al_at 


-1.725898 


GPM6B 


Glycoprotein M6B 


NM_001001994 


Ssc.28842.1.Al_at 


-1.685047 


CCDC113 


Coiled-coil domain containing 113 


NM_014157 


Ssc.7917.1.Al_at 


-1.680845 


Null 


MI-P-E6-ahn-c-09- 1 -UM.s 1 M1-P-E6 Sus scrofa cDNA clone 


BQ599768.1 


Ssc.29910.1.Al_at 


-1.633574 


MTMR2 


Myotubularin related protein 2 


NM_016156 


Ssc.6616.1.Sl_at 


-1.629585 


Clor/92 


Chromosome 1 open reading frame 92 


NM_ 144702 


Ssc.28950.1.Al_at 


-1.627009 


EX0C6 


Exocyst complex component 6 


NM_019053.4 
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Table 5. Functional annotation clusters of up-regulated genes in uterine endometrium on D12 of pregnancy compared with those on D12 
of the estrous cycle 



Annotation cluster^ 


Enrichment score'' 


Biological terms' 


1 


2.56 


Membrane fraction (12), cell fraction (14), insoluble fraction (12) 


2 


1.04 


Endoplasmic reticulum part (6), endoplasmic reticulum membrane (5), endoplasmic reticulum (10), nuclear 
envelope-endoplasmic reticulum network (5), microsome (4), vesicular fraction (4), organelle membrane 
(9), endomembrane system (6) 


3 


0.88 


Blood vessel morphogenesis (4), blood vessel development (4), vasculature development (4) 


4 


0.86 


Cellular macromolecule catabolic process (9), proteolysis (11), macromolecule catabolic process (9), 
proteolysis involved in cellular protein catabolic process (7), cellular protein catabolic process (7), protein 
catabolic process (7), ubiquitin- protein ligase activity (3), small conjugating protein ligase activity (3), 
Ugase activity, forming carbon-nitrogen bonds (3), ubiquitin-dependent protein catabolic process (3), 
modification-dependent macromolecule cataboUc process (6), acid- amino acid ligase activity (3), 
modification-dependent protein catabolic process (6) 


5 


0.80 


Intracellular organelle lumen (14), organelle lumen (14), nuclear lumen (10), membrane-enclosed lumen 

(14), nucleoplasm (8) 


6 


0.72 


Proteolysis (11), metalloendopeptidase activity (3), endopeptidase activity (5), negative regulation of 
apoptosis (5), negative regulation of programmed cell death(5), negative regulation of cell death (5), 
metallopeptidase activity (3), peptidase activity, acting on L-amino acid peptides (5), peptidase activity (5), 
regulation of apoptosis (6), regulation of programmed cell death (6), regulation of cell death (6) 


7 


0.72 


Phosphatase activity (5), phosphoprotein phosphatase activity (4), phosphate metabolic process (9), 
dephosphorylation (4), phosphorus metabolic process (9), protein amino acid dephosphorylation (3), protein 
serine/threonine kinase activity(3), phosphorylation (5), protein amino acid phosphorylation (5), protein 
kinase activity (3) 


8 


0.71 


Positive regulation of multicellular organismal process (4), positive regulation of cell proliferation (5), 
regulation of growth (4) 


9 


0.65 


Cell cycle (8), mitotic cell cycle (5), cell cycle process (5), cell cycle phase (4) 


10 


0.65 


Glycosaminoglycan binding (3), polysaccharide binding (3), pattern binding (3), carbohydrate binding (3) 



" Top ten of amotation clusters identified by official gene symbols of up-regulated gene list on D12 of pregnancy compared with those on D12 of the estrous cycle through 
DAVID. 



^ Enrichment score represents significance of each cluster annotation and the relatedness of the terms and the genes associated with terms. Enrichment score is calculated by 

overall EASE scores (the modified Fisher extract p value) of each term member. 
^ The gene ontology (GO) terms are group of terms having similar biological meaning with respect to biological process, cellular component, and molecular function. GO 

functional annotation term (FAT) was used for the analysis because GO term cannot univereally define the specificity of a given term. The number in parenthesis indicate the 

number of differentially expressed genes contribute to the clustered term. 



Table 6. Functional annotation clusters of down-regulated genes in uterine endometrium on D12 of pregnancy compared with those on 

D12 of the estrous cycle 

Annotation cluster' Enrichment score'' Biological tcinis' 

1 2.04 Plasma membrane part (22), integral to plasina inembrane (11), intrinsic to plasma membrane (11) 

2 2.03 Cell adhesion (11), biological adhesion (11), cell-cell adhesion (4) 

3 1.50 Leukocyte activation (7), cell activation (7), lymphocj'te activation (6), T cell activation (4), immune 

effector process (4), regulation of leukocyte activation (4), regulation of cell activation (4), leukocyte 
mediated immunity (3), cell surface (5), positive regulation of response to stimulus (4), positive regulation 
of immune system process (4), immune response (7), external side of plasma membrane (3) 

4 1.48 Immune system development (7), hemopoiesis (6), hemopoietic or lymphoid organ development (6), 

myeloid cell differentiation (3), protein homodimerization activity (4), nucleoplasm part (3) 

5 1.46 Cell junction (7), adherens junction (4), anchoring junction 4) 

6 1.45 Biopolymer glycosylation (4), glycoprotein biosynthetic process (4), glycoprotein metaboUc process (4), 

glycosylation (4), protein amino acid glycosylation (4) 

7 0.98 Apoptosis (7), programmed cell death (7), cell death (7), death (7), negative regulation of cell proliferation 

(4) 

8 0.83 DNA metabolic process (8), response to DNA damage stimulus (5), DNA repair (4), cellular response to 

stress (5), DNA binding (7) 

9 0.83 Induction of apoptosis (5), induction of programmed cell death (5), positive regulation of apoptosis (5), 

positive regulation of programmed cell death (5), positive regulation of cell death (5), regulation of 
apoptosis (7), regulation of prograimned cell death (7), regulation of cell death (7) 
10 0.76 Response to ionizing radiation (3), response to radiation (3), response to abiotic stimulus (3) 

" Top ten of annotation clusters identified by official gene symbols of up-regulated gene list on D15 of pregnancy compared with those on D15 of the estrous cycle. 
''Enrichment score represents significance of each cluster annotation and the relatedness of the terms and the genes associated with terms. Enrichment score is calculated by 
overall EASE scores (the modified Fisher extract p value) of each term member. 

The gene ontology (GO) terms are group of terms having similar biological meaning with respect to biological process, cellular component, and molecular function. GO 
functional annotation term (FAT) was used for the analysis because GO term cannot universally define the specificity of a given term. The number in parenthesis indicate the 
number of differentially expressed genes contribute to the clustered term. 
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Table 7. Functional annotation clusters of up-regulated genes in uterine endometrium on D15 of pregnancy compared with those on D15 
of the estrous cycle 



Annotation 








Enrichment score'' 


Biological terms'^ 


1 


3.46 


Nuclear speck (7), nuclear body (8), RNA splicing (9) 


2 


2.21 


Nucleotide binding (30), purine ribonucleotide binding (25), ATP binding (21), ribonucleotide binding (25), purine 






nucleotide binding (25), adenyl ribonucleotide binding (21), adenyl nucleotide binding (21), purine nucleoside binding 






(21), nucleoside binding (21) 


3 


2.12 


Nuclear body (8), nucleoplasm part (12), nucleoplasm (15), membrane-enclosed lumen (24), nuclear lumen (20), 






intracellular organelle lumen (23), organelle lumen (23), nucleolus (9) 


4 


1.51 


RNA splicing (9), mrna metabolic process (9), RNA processing (11), mrna processing (8), spHceosome assembly (3), RNA 






splicing, via transesterilication reactions (5), RNA splicing, via transesterification reactions with bulged adenosine as 






nucleophile (5), nuclear mrna splicing, via spliceosome (5), RNA binding (11), ribonucleoprotein complex assembly (3), 






ribonucleoprotein complex biogenesis (4), cellular macro molecular complex assembly (5), cellular macromolecular 






complex subunit organization (5) 


5 


1.15 


Nuclear body (8), nuclear export (4), nucleocytoplasmlc transport (4), nuclear transport (4), mrna transport (3), 






establishment of RNA localization (3), nucleic acid transport (3), RNA transport (3), RNA localization (3), nucleobase. 






nucleoside, nucleotide and nucleic acid transport (3), intracellular transport (8) 


6 


1.06 


Receptor-mediated endocytosis (4), membrane organization (8), membrane invagination (4), endocytosis (4), vesicle- 






mediated transport (6) 


7 


0.96 


Helicase acti\'ity (6). ATP-dependent helicase activity (3), purine NTP-dependent helicase activity (3), atpase activity. 






coupled (4), atpase activity (4) 


8 


0.95 


Enzyme binding (8), protein kinase binding (4). kinase binding(4) 


9 


0.94 


Macromolecular complex assembly (11), macromolecular complex subunit organization (11), protein complex assembly 






(8), protein complex biogenesis (8), cellular macromolecular complex assembly (5), cellular macromolecular complex 






subunit organization (5), protein oligomerization (3) 


10 


0.75 


Protein serine/threonine kinase activity (7), enzyme linked receptor protein signaling pathway (6), phosphate metabolic 






process (11), phosphorus metabolic process (11), protein kinase activity (8), protein amino acid phosphorylation (8), 






phosphorylation (9) 


Top ten of annotation ckistcrs identified by official 


gene symbols of up-regulated gene list on D12 of pregnancy compared with those on D12 of the estrous cycle. 


Enrichment score represents significance of each cluster annotation and the relatedness of the terms and the genes associated with terms. Enrichment score is calculated by overall EASE 


scores (the modified Fisher extract p value) of each term member. 


" The gene ontology (GO) terms are group of terms having similar biological meaning with respect tt) biological process, cellular component, and molecular function. GO functional 


annotation term (FAT) was used for the analysis because GO term cannot universally define the specificity of a given term. The number in parenthesis indicate the number of differentially 


expressed genes contribute to the clustered term. 




Table 8. Functional annotation clusters of down-regulated genes in uterine endometrium on D 15 of pregnancy compared with those on 


D15 of the estrous cycle 




Annotation 






cluster^ 


Enrichment score'' 


Biological terms'^ 


1 


1.68 


Adenyl ribonucleotide binding (15), purine ribonucleotide binding (17), purine nucleotide binding (17), ribonucleotide 






binding (17), nucleoside binding (15), adenyl nucleotide binding (15), purine nucleoside binding (15), ATP binding (14), 






nucleotide binding (18) 


2 


1.50 


Cytoskeletal part (11), microtubule cytoskeleton (7), intracellular non-membrane- bounded organelle (18), non-membrane- 






bounded organelle (18), centrosome (4), cytoskeleton (13), microtubule organizing center (4) 


3 


1.50 


Cytoskeletal part (11), dynein complex (3), microtubule associated complex (4), microtubule-based movement (4), 






microtubule cytoskeleton (7), motor activity (3), microtubule motor activity (3), microtubule (4), microtubule-based 






process (4) 


4 


1.20 


Cytoskeletal part (11), regulation of actin filament length (3), regulation of cellular component size (5), regulation of actin 






cytoskeleton organization (3), regulation of actin filament-based process (3), regulation of organelle organization (4), actin 






cytoskeleton (4), cytoskeletal protein binding (6), regulation of cytoskeleton organization (3), actin binding (4) 


5 


1.17 


Regulation of neuron differentiation (4), regulation of neurogenesis (4), regulation of nervous system development (4), 






regulation of cell development (4), neuron differentiation (4) 


6 


0.92 


Protein amino acid dephosphorylation (4), dephosphorylation (4), phosphorus metabolic process (9), phosphoprotein 






phosphatase activity (4), phosphate metabolic process (9), protein tyrosine phosphatase activity (3), phosphatase activity 






(4), protein kinase activity (5), protein serine/threonine kinase activity (4), protein amino acid phosphorylation (5), 






phosphorylation (5) 


7 


0.84 


Central nervous system neuron axonogenesis (3), central nervous system neuron development (3), central nervous system 






neuron differentiation (3), cell projection organization (5), cell motion (5), axonogenesis (3), cell morphogenesis involved 






in neuron differentiation (3), neuron projection morphogenesis (3), cell projection morphogenesis (3), cell morphogenesis 






involved in differentiation (3), cell part morphogenesis (3), neuron differentiation (4), neuron projection development (3), 






cell morphogenesis (3), cellular component morphogenesis (3), neuron development (3) 


8 


0.82 


Regulation of synaptic transmission (3), regulation of transmission of nerve impulse (3), regulation of neurological system 






process (3), regulation of system process (4) 


9 


0.75 


Cell projection organization (5), cytoskeletal protein binding (6), cytoskeleton organization (3) 


10 


0.71 


Learning or memory (3), neurological system process (10), behavior (5), cognition (4) 



* Top ten of annotation clusters identified by official gene symbols of up-regulated gene hst on D15 of pregnancy compared with those on D15 of the estrous cycle. 

''Enrichment score represents significance of each cluster annotation and the relatedness of the terms and the genes associated with terms. Enrichment score is calculated by overall EASE 
scores (the modified Fisher extract p value) of each term member. 



The gene ontology (GO) terms are group of terms having similar biological meaning with respect to biological process, cellular component, and molecular function. GO functional 
annotation term (FAT) was used for the analysis because GO term cannot universally define the specificity of a given term. The number in parenthesis indicate the number of differentially 
expressed genes contribute to the clustered term. 



